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HIGHLIGHTS 


•  We  modeled  the  electrical  and  thermal  behavior  of  the  Li-ion  battery. 

•  We  analyzed  the  dynamic  behavior  of  the  cell  using  SFUDS  cycle. 

•  We  carried  out  experimental  study  to  validate  the  simulation  results. 

•  We  studied  the  thermal  response  of  the  battery  pack  under  different  driving  cycle. 

•  Heat  generation  of  the  battery  pack  is  the  highest  for  US06  cycle. 
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Lithium  ion  batteries  offer  an  attractive  solution  for  powering  electric  vehicles  due  to  their  relatively  high 
specific  energy  and  specific  power,  however,  the  temperature  of  the  batteries  greatly  affects  their  per¬ 
formance  as  well  as  cycle  life.  In  this  work,  an  empirical  equation  characterizing  the  battery’s  electrical 
behavior  is  coupled  with  a  lumped  thermal  model  to  analyze  the  electrical  and  thermal  behavior  of  the 
18650  Lithium  Iron  Phosphate  cell.  Under  constant  current  discharging  mode,  the  cell  temperature  in¬ 
creases  with  increasing  charge/discharge  rates.  The  dynamic  behavior  of  the  battery  is  also  analyzed 
under  a  Simplified  Federal  Urban  Driving  Schedule  and  it  is  found  that  heat  generated  from  the  battery 
during  this  cycle  is  negligible.  Simulation  results  are  validated  with  experimental  data.  The  validated 
single  cell  model  is  then  extended  to  study  the  dynamic  behavior  of  an  electric  vehicle  battery  pack.  The 
modeling  results  predict  that  more  heat  is  generated  on  an  aggressive  US06  driving  cycle  as  compared  to 
UDDS  and  HWFET  cycle.  An  extensive  thermal  management  system  is  needed  for  the  electric  vehicle 
battery  pack  especially  during  aggressive  driving  conditions  to  ensure  that  the  cells  are  maintained 
within  the  desirable  operating  limits  and  temperature  uniformity  is  achieved  between  the  cells. 

©  2013  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  world  relies  heavily  on  fossil  fuel  to  meet  the  daily  power 
demands,  ranging  from  electricity  generation  to  transportation.  In 
2009,  the  logistics  sector  had  contributed  to  61.7%  of  the  total  world 
oil  consumption  and  23%  of  the  total  world  CO2  emission  [1].  With 
advances  in  battery  technology,  concerns  on  global  warming  and 
increasing  fuel  prices,  automotive  manufacturers  are  forced  to  shift 
their  attention  to  electric  vehicles  (EVs).  EVs  are  more  energy 
efficient  and  cleaner  than  conventional  Internal  Combustion 
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Engine  (ICE)  vehicles  and  are  projected  as  the  most  sustainable 
solutions  for  the  future  transport  [1-3],  however,  the  success  of  EVs 
depends  on  the  development  of  the  battery. 

Lithium-ion  batteries  provide  an  attractive  solution  for  EVs  due 
to  its  high  power  and  energy  density,  however,  thermal  issues  in  Li- 
ion  batteries  have  to  be  addressed  to  make  them  safer,  reliable  and 
last  longer  for  high  power  applications.  Bandhaeur  et  al.  4]  have 
provided  a  detailed  review  of  the  thermal  issues  in  Li-ion  batteries. 
Safety,  cycle  life  and  capacity  retention  are  some  of  the  major  as¬ 
pects  affected  by  the  operating  temperature.  Hence,  a  control  and 
management  of  the  thermal  envelope  is  required  to  have  a  safer 
and  reliable  operation  of  a  Li-ion  battery. 

Various  test  procedures  used  to  investigate  the  performance, 
reliability  and  safety  of  Li-ion  batteries  are  documented  in  ISO 
12405-1/2,  IEC  62660-1/2  [3].  In  order  to  design  a  realistic  thermal 
management  system  for  EVs,  it  is  also  important  to  characterize  the 


232 


L.H.  Saw  et  al.  /  Journal  of  Power  Sources  249  (2014)  231-238 


Nomenclature 

Rc 

terminal  contact  resistance,  0 

Tdis 

discharge  time,  h 

A 

exponential  voltage,  V 

^module 

module  temperature,  K 

As 

external  surface  area  of  the  cell,  m2 

Tsurf 

surface  temperature  of  battery,  K 

71  ITT  f 

cross  section  flow  area  for  cooling  air  per  module,  m2 

Toe 

free  stream  temperature  of  air,  K 

71ms 

total  module  surface  area  exposed  to  cooling  air,  m2 

t 

time,  s 

B 

exponential  zone  time  constant,  (Ah)-1 

VW 

voltage,  V 

cP 

specific  heat  capacity  of  the  battery,  J  kg-1  K-1 

P 

density  of  the  battery,  kg  m-3 

Cbp 

theoretical  capacity  of  the  battery,  Ah 

Pair 

density  of  the  air,  kg  m-3 

E 

emissivity 

^sb 

Stefan— Boltzmann  constant,  W  m-2  K-4 

Eo 

constant  voltage,  V 

h 

convective  heat  transfer  coefficient,  W  m-2  K-1 

Subscript 

fdis 

discharge  current,  A 

exp 

exponential 

i 

current,  A 

nom 

nominal 

i* 

filtered  current,  A 

full 

fully  charged 

I< 

polarization  constant,  V(Ah)-1  or  polarization 

max 

maximum 

resistance,  Q 

m 

mass  flow  rate  of  air,  kg  s-1 

Superscripts 

a 

capacity,  Ah 

k 

Peukert  constant 

R 

internal  resistance,  Q 

thermal  phenomena  of  Li-ion  cell  for  the  required  transient  power 
response.  Among  the  various  driving  cycles,  Simplified  Federal 
Urban  Driving  Schedule  (SFUDS)  has  been  widely  used  to  evaluate 
the  power  delivery  capability  and  cycle  life  of  the  cell  at  laboratory 
level  [5-7],  however,  the  facility  needed  to  carry  out  the  testing  of 
the  EV  battery  packs  such  as  high  power  programmable  battery 
cycler  is  always  expensive  and  requires  several  hundred  hours  of 
testing.  The  battery  pack  may  not  be  comprehensively  tested  due  to 
the  limitations  of  the  battery  cycler. 

Numerical  modeling  could  be  used  to  overcome  the  limita¬ 
tions  in  battery  testing.  Numerical  modeling  not  only  helps  to 
improve  the  understanding  of  the  battery  operating  mechanism 
but  also  provides  internal  information  that  are  difficult  to  obtain 
through  experiments  such  as  electrochemical  reaction  rates  in¬ 
side  the  cell,  heat  generation,  temperature  distribution,  voltage 
distribution,  current  distribution,  etc.  Various  mathematical 
models  have  been  used  to  investigate  the  thermal  response  of  the 
battery  such  as  empirical  equations  [7],  electrochemical  models 
[8-12],  RC  models  [13-15],  and  lumped  parameter  models 
[16,17].  Although  electrochemical  models  can  predict  the  aging 
and  thermal  behavior  of  the  Li-ion  battery,  coupled  time  variant 
spatial  partial  differential  equations  make  them  complex  and 
their  solution  demands  extensive  computational  resources  18]. 
Besides,  most  of  the  studies  only  presented  a  numerical  simula¬ 
tion  results  and  doesn’t  compare  with  the  experimental  work 
[4,19-22].  Various  types  of  thermal  management  systems  such  as 
air  cooling  [23-25],  liquid  cooling  [26,27]  and  phase  change 
material  [28-30]  technologies  have  been  proposed  to  maintain 
the  batteries  at  the  optimum  operating  temperature  at  both  cell 
and  pack  levels  [31]. 

In  view  of  the  challenges  in  experimental  testing  and  detailed 
modeling,  the  objective  of  this  work  is  twofold.  First,  an  empirical 
equation  coupled  with  a  lumped  thermal  model  has  been  used  to 
predict  the  cell  voltage,  heat  generation,  temperature  rise  of  the  cell 
during  constant-current  discharging  and  SFUDS  cycle  for  an  18650 
Lithium  Iron  Phosphate  (LFP)  cell  and  is  validated  with  experi¬ 
ments;  and  second,  to  apply  the  validated  single  cell  model  to 
investigate  the  thermal  response  of  the  battery  pack  of  a  converted 
EV  under  Urban  Dynamometer  Driving  Schedule  (UDDS),  Flighway 
Fuel  Economy  Driving  Schedule  (FIWFET)  and  US06  Supplemental 
Federal  Test  Procedure  (SFTP)  driving  cycles.  The  results  are  dis¬ 
cussed  in  terms  of  the  total  heat  generated  during  these  driving 


cycles  and  the  evolution  of  the  battery  pack  temperature  under  a 
forced  convection  cooling  system. 

2.  Mathematical  model 

2.1.  The  battery  model 

A  battery  model  is  needed  to  define  its  voltage  in  terms  of 
current  and  state  of  charge  (SOC).  In  this  study,  modified  Shepherd 
model  has  been  employed  to  represent  the  voltage  dynamics  of  the 
LFP  cell  [32-34  .  A  typical  discharge  curve  of  the  Li-ion  battery  is 
shown  in  Fig.  2.  The  discharge  curve  of  the  Li-ion  battery  can  be 
divided  into  three  sections.  The  first  section  represents  the  expo¬ 
nential  potential  drop  of  the  cell  during  initial  discharge.  The  sec¬ 
ond  section  represents  the  amount  of  charge  that  can  be  extracted 
from  the  cell  before  reaching  the  nominal  voltage  of  the  cell.  The 
last  section  represents  the  total  discharge  of  the  cell,  when  the 
voltage  of  the  cell  drops  rapidly  to  the  cut  off  voltage.  The  modified 
Shepherd  equation  for  charging  and  discharging  is  given  by  Eqs.  (1 ) 
and  (2)  respectively  [32-34  .  It  is  assumed  that  the  internal  resis¬ 
tance  of  the  cell  is  constant  throughout  the  charging  and  dis¬ 
charging  cycle  and  doesn’t  change  with  the  /t-rates.  The 
temperature  effect  on  the  battery  model  behavior  is  neglected  and 
the  model  parameters  for  discharging  and  charging  are  identical. 
Charging  (i*  <  0) 

VW  =  E0-Ri-K  .t_g1Q  ■  i*  -  K  it  +  Aexp(  -  B  •  it) 

(1) 

Discharging  (i*  >  0) 

VW  =  E0-  R  i  ~  I<q% Tt-f  -l<((^T_yt  +  Aexp(-B  it) 

(2) 

The  voltage  of  the  cell  at  a  fully  charged  state  is  defined  in  Eq.  (3) 

[33,34], 

Viuu  =  E0-Ri  +  A  (3) 

The  voltage  at  the  exponential  section  is  defined  in  Eq.  (4) 

[33,34], 
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Fig.  1.  Testing  cycle  for  single  cell  and  EV  battery  pack. 


Vexp  —  Eq  —  K  - (Q.  exp  +  i)  —  R  i  Aexp  (  — - Q.  exp  ) 

—  kL  exp  yCi  exp  J 

(4) 

The  voltage  of  the  cell  at  the  nominal  zone  is  defined  in  Eq.  (5) 

[33,34]. 


^nom  —  E0-I<  — - — - (Q.  nom  +  0  —  Ri 

Cl  C:  nom 

+  Aexp  (  — - Q.nom  )  (5) 

\\Lexp  J 


Fig.  2.  Typical  discharge  characteristic  of  Li-ion  battery  [33,34]. 


2.2.  Thermal  model 

It  is  assumed  that  the  resistance  to  conduction  within  the  bat¬ 
tery  is  much  less  than  the  resistance  to  convection  outside  the 
battery  surface.  Therefore,  the  temperature  inside  the  cell  is  fairly 
uniform  [35  .  Besides,  it  is  also  assumed  that  the  physical  properties 
of  the  cell  are  uniform  and  are  not  affected  significantly  by  tem¬ 
perature.  The  battery  model  is  coupled  with  a  lumped  thermal 
model  to  determine  the  heat  generation  and  the  average  surface 
temperature  of  the  cell.  A  general  energy  balance  equation  used  to 
model  the  battery  system  is  defined  in  Eq.  (6)  [7]. 

pCpT  =  Q  sen  “  hAs  (T*urf  -  )  -  Easb  (7?urf  -It)  (6) 

The  heat  generation  in  the  Li-ion  battery  constitutes  of  ohmic 
heat,  irreversible  heat  and  reversible  heat  [36].  Reaction  heat  or 
irreversible  heat  generation  is  due  to  the  transfer  of  electrons  to  or 
from  the  electrode  during  the  electrochemical  reaction,  reversible 
heat  generation  is  because  of  the  entropy  changes  of  cathode  and 
anode,  and  ohmic  heat  generation  is  due  to  the  ohmic  resistance  of 
the  solid  active  materials  and  electrolyte.  The  entropy  change  in  the 
electrodes  is  related  to  the  change  of  their  equilibrium  potential 
with  temperature  (dH/dT)  and  this  varies  with  SOC.  Here,  this 
relation  is  adopted  from  the  work  of  Forgez  et  al.  [37].  Heat  gen¬ 
eration  due  to  contact  resistance  is  also  added  to  the  overall  heat 
generation  within  the  cell  as  given  in  Eq.  (7)  [38]. 

Qgen  =  ‘'(Eel  -  V^tt  +  +  fRc  (7) 

2.3.  Battery  pack  thermal  model 

Different  types  of  cooling  systems  will  influence  the  perfor¬ 
mance  and  cost  of  the  battery  pack  thermal  management  system. 
The  heat  transfer  medium  could  be  air,  liquid,  heat  pipes  and  phase 
change  material  (PCM),  or  a  combination  of  them.  The  architecture 
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Exhaust  Air  Cross  Section 

Outlet :  Tair  for  internal  flow 


Fig.  3.  Schematic  of  active  air  cooling  system  of  an  EV  battery  pack  [13]. 


of  cooling  strategy  includes  series  flow,  parallel  flow  and  combi¬ 
nation  of  series  and  parallel  flow.  The  selection  of  the  cooling  sys¬ 
tem  depends  on  the  constraints  of  the  vehicle,  installation  costs  and 
the  external  environment.  In  the  extreme  environment  and  work¬ 
ing  under  heavy  duty  cycles,  an  active  cooling  system  is  preferred  to 
offer  more  effective  thermal  management.  In  this  study,  parallel  air 
flow  strategy  is  adopted  to  investigate  the  thermal  response  of  the 
battery  pack  with  a  certain  number  of  modules.  Each  module  will 
have  the  same  inlet  air  temperature  that  will  result  in  more  uni¬ 
form  pack  temperature  [13  and  hence  it  is  sufficient  to  study  a 
single  module  in  the  entire  pack.  A  schematic  diagram  of  the  bat¬ 
tery  pack  and  the  air  flow  is  shown  in  Fig.  3.  The  heat  transfer  co¬ 
efficient  for  natural  convection  and  forced  convection  is  defined  in 
Eq.  (8)  [13,39]. 

h  =  |  ^ forced  =  30  (  '  rmodule  >  35  "C  (g) 

'  ^natural  m  4  ^module  —  35  C 


The  surface  temperature  of  the  battery  was  measured  using 
twelve  thermocouples  (T-type)  attached  to  different  locations  on 
the  cells.  Three  thermocouples  were  attached  in  the  axial  direction 
and  four  sides  of  the  battery  surface  as  shown  in  Fig.  4.  Measure¬ 
ment  of  battery  surface  temperature  during  charging  and  dis¬ 
charging  at  different  /t-rates  was  done  at  room  temperature  of  25  °C 
under  natural  convection.  All  the  tests  were  repeated  three  times 
and  the  average  value  was  taken.  Temperature  readings  were 
recorded  using  the  HP  34970A  data  acquisition  system.  The  specific 
heat  capacity  of  the  cell  is  measured  using  adiabatic  accelerating 
rate  calorimeter  (THT  ARC).  Vehicle  specific  parameters  used  for 
driving  cycle  simulations  are  tabulated  in  Table  2. 

2.5.  Numerical  procedures 

The  battery  and  thermal  model  equations  are  solved  simulta¬ 
neously  in  the  Matlab-Simulink  2011b.  Parameters  extracted  in 
Section  2.4  are  entered  into  the  battery  model  in  Simulink  to 
investigate  the  electrical  and  thermal  response  of  the  cell  during 
constant  current  discharging  and  SFUDS.  Outputs  of  the  simulation 
model  are  voltage,  current,  heat  generation  and  surface  tempera¬ 
ture  of  the  cell.  Heat  is  dissipated  from  the  cell  by  natural  con¬ 
vection  and  radiation.  The  model  was  validated  using  constant 
current  discharging  experiments  ranging  from  1  Jt  (1.3  A)  to  3  Jt 
(3.9  A)  and  SFUDS  dynamic  power  profile.  Under  SFUDS  test,  the 
cell  is  subjected  to  continuous  charging  and  discharging  until  its 
SOC  reaches  10%.  This  validated  model  was  then  used  to  study  the 
thermal  response  of  the  battery  pack  for  a  converted  EV  using 
Hyundai  Trajet  under  UDDS,  HWFET  and  US06  driving  cycles  that 
are  illustrated  in  Fig.  1. 


2.4.  Experimental  setup  and  parameters  extraction 

Commercial  1.3  Ah  18650  cells  with  graphite  anode  coated  on 
the  copper  current  collector,  LFP  cathode  coated  on  the 
aluminum  current  collector,  filled  with  an  electrolyte  of  LiPF6  in 
EC:DEC  1:1  and  Polyvinylidene  Fluoride  (PVDF)  separator  were 
used  for  experiments.  The  charging  and  discharging  of  the  bat¬ 
tery  was  conducted  using  a  battery  cycler  (Maccor  Instrument 
4000).  According  to  the  procedure  followed  to  Tremblay  et  al. 
[33,34],  the  cell  is  discharged  at  0.2  Jt  and  the  parameters  Fo,  K  A 
and  B  are  extracted  from  the  discharge  curve  shown  in  Fig.  5(a). 
The  cut  off  voltage  for  constant  current  discharging  was  2.3  V. 
Charging  was  done  in  two  modes:  constant  current  charging  until 
4.2  V  followed  by  constant  voltage  charging  until  the  current 
dropped  to  0.1  A.  Before  discharging,  the  cell  is  charged  at  0.1  Jt 
and  followed  by  1  h  of  rest.  Slow  rate  of  charging  is  necessary  to 
ensure  that  the  chemical  process  within  the  cell  occurs  at  the 
similar  rate  to  the  transfer  of  electric  energy.  Internal  resistance 
(R)  of  the  battery  is  measured  using  an  impedance  analyzer 
(Solartron  analytical  1400).  The  values  of  the  extracted  parame¬ 
ters  and  other  dimensions  used  in  the  simulation  are  provided  in 
Table  1. 


Fig.  4.  Experimental  setup  for  temperature  measurement  and  the  location  of  ther¬ 
mocouples  attached  to  the  surface  of  the  cell. 


Time,  min 

Fig.  5.  (a)  Model  calibration:  Cell  voltage  during  discharge  at  0.2  It  to  extract  the  model 
parameters,  (b)  Validation  of  battery  model:  Cell  voltage  during  discharge  at  1,  2,  and  3 
It- rates. 
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Table  1 

LFP  cell  parameters. 


Parameter 

Value 

Parameter 

Value 

Nominal  voltage,  V 

3.0 

Eo,  V 

3.21 

Nominal  capacity,  Ah 

1.3 

R,  Q 

0.03 

Weight,  kg 

0.040 

I<,  Q  or  V(Ah)-1 

0.0119 

Density,  kg  m-3  [40] 

2700 

A,  V 

0.2711 

Cathode  material 

LiFeP04 

B,  (Ah)"1 

152.130 

Anode  material 

Graphite 

Reference  temperature, 

Tref(K) 

298.15 

Cell  length,  m 

0.065 

Cell  diameter,  m 

0.018 

Specific  heat  capacity, 

J  kg"1 K-1 

900 

Emissivity  of  heat 
shrink  wrapping,  e  [41  ] 

0.95 

Stefan— Boltzmann 
constant,  a  (W  m 
-2  K4) 

5.67  x  10“8 

Heat  transfer 
coefficient  for  constant 
current  discharging, 

W  ur2  K1  (estimated) 

8 

3.  Result  and  discussion 

3.1.  Validation  of  the  cell  potential 

Discharge  characteristics  of  the  cell  predicted  by  the  battery 
model  and  experimental  data  are  provided  in  Fig.  5(a).  The  average 
squared  error  of  the  simulated  and  measured  voltage  is  about 
0.000184.  This  is  followed  by  comparing  the  predicted  results  of  the 
model  with  the  experimental  discharge  curves  at  1, 2,  and  3  Jt- rates 
as  shown  in  Fig.  5(b).  The  model  predictions  match  well  with  the 
experiments  with  the  average  squared  error  of  0.000482,  0.000184 
and  0.000752  respectively  for  the  various  discharge  rates.  At  3  Jt  of 
constant  current  discharging,  the  maximum  capacity  of  the  cell  was 
reduced  by  8%  less  than  the  capacity  at  0.2  Jt  (1.3  Ah).  The  reduction 
in  the  cell  capacity  at  high  Jrrates  can  be  explained  by  Peukert 
equation  as  in  Eq.  (9).  The  discharge  capacity  of  the  cell  reduces  as 
the  discharge  current  increases  [3]: 

Cbp  =  Tdis-4  (9) 

3.2.  Evolution  of  the  cell  temperature  and  heat  generation 

The  average  measured  and  simulated  temperature  rise  of  the 
cell  surface  at  different  Jt-rates  of  discharging  is  shown  in  Fig.  6.  The 
temperature  of  the  cell  has  a  positive  relationship  with  the  /t-rates. 
The  sharp  rise  of  the  temperature  towards  the  end  of  discharge  was 
probably  due  to  the  polarization  effect  of  the  cell.  The  simulated 
results  agree  well  with  the  experimental  results.  The  average 
squared  error  between  the  simulation  results  and  the  experiments 
at  1,  2,  and  3  Jt-rates  of  discharging  are  0.00961,  0.00976  and 
0.00831  respectively.  The  measured  average  temperature  rise  of  the 
cell  at  the  end  of  discharging  at  1,  2,  and  3  Jt-rates  is  6.4  °C,  13.9  °C 
and  20.0  °C  respectively.  At  low  Jt-rates  of  discharging,  heat 
generated  from  the  cell  can  be  effectively  dissipated  by  natural 
convection  and  only  a  temperature  rise  of  less  than  10  °C  is 
observed.  On  the  other  hand,  large  amount  of  heat  was  generated  at 
3  It  and  natural  convection  is  not  sufficient  to  dissipate  it  and  keep 
the  battery  within  the  recommended  operating  temperature  range. 
The  total  heat  generation  predicted  by  modeling  at  different  Jt- rates 
during  discharging  is  shown  in  Fig.  7.  Large  amount  of  heat  is 
generated  towards  the  end  of  the  discharge.  The  maximum  amount 
of  heat  generated  for  1,  2,  and  3  /t-rates,  is  0.59  W,  1.47  W  and 
3.14  W  per  cell  respectively. 

3.3.  Dynamic  behavior  of  the  cell  under  SFUDS 

In  SFUDS  testing  cycle,  the  battery  is  operated  under  a  repeated 
360  s  cycle  of  charging  and  discharging  at  certain  specific  power  as 


shown  in  Fig.  1.  Variation  of  the  current  and  voltage  during  the  first 
cycle  of  the  SFUDS  profile  is  shown  in  Fig.  8(a)  and  (b)  respectively. 
Maximum  current  for  discharging  is  about  0.97  A  corresponding  to 
a  specific  power  of  79  W  kg-1.  The  average  squared  error  in  the 
experimental  data  and  simulated  results  of  current  is  0.000062. 
The  over  potential  of  the  cell  during  the  first  cycle  is  high  ( ~3.3  V) 
and  amount  of  energy  being  discharged  from  the  cell  is  about  1%. 
The  average  squared  error  between  the  experimental  data  and 
simulated  results  of  over  potential  is  found  to  be  0.0000046. 

Thermal  response  for  the  battery  is  studied  for  more  number  of 
repeated  cycles  until  it  reaches  the  cut-off  voltage  and  the  pre¬ 
dicted  average  heat  generation  data  is  provided  in  Fig.  9.  The 
average  heat  generation  rate  of  the  cell  during  600  min  of  the 
SFUDS  cycle  is  relatively  small  as  compared  to  the  constant  current 
discharging  and  can  be  neglected.  The  average  heat  generated  from 
the  first  cycle  of  SFUDS  is  about  -4.6  x  10  4  mj.  Endothermic  heat 
generation  is  observed  during  the  initial  75  min  of  the  cycle  and 
heat  generation  rate  is  gradually  reduced  as  the  dynamic  charging 
and  discharging  process  proceeds.  Towards  the  end  of  the  SFUDS 
cycle,  heat  generated  from  the  cell  contributes  to  the  increase  of  the 
overall  cell  temperature. 

The  comparison  of  experimental  data  and  model  predictions  for 
the  voltage  and  temperature  profile  of  the  cell  under  600  min  of 
SFUDS  cycle  is  shown  in  Fig.  10(a)  and  (b)  respectively.  In  general, 
the  voltage  of  the  cell  decreases  with  the  number  of  cycles  and 
oscillate  in  phase.  As  shown  in  Fig.  10(a),  there  is  only  a  small  dif¬ 
ference  between  the  measured  and  simulated  cell  voltages.  Sub¬ 
stantial  amount  of  heat  is  generated  towards  the  end  of  the  SFUDS 
cycle  as  compared  to  the  initial  6  min  of  the  cycle.  Therefore,  the 
temperature  of  the  cell  slowly  increases  at  the  end  of  the  cycle  as 
shown  in  Fig.  10(b).  The  heat  generated  in  the  cell  can  be  effectively 
dissipated  by  natural  convection  and  cell  is  kept  within  safer 
operating  limits.  This  validated  model  can  now  be  extended  to  a 
battery  pack  and  employed  to  study  the  thermal  response  of  it 
under  various  driving  cycles  as  discussed  in  the  next  section. 

3.4.  Thermal  responses  of  the  battery  pack  under  different  driving 
cycle 

The  thermal  response  of  the  battery  pack  designed  for  the 
converted  EV  under  different  driving  cycles  is  provided  in  Figs.  11- 
13.  A  thermostat  feature  is  also  included  in  the  Simulink  model 
with  a  set  temperature  of  35  °C  i.e.,  air  flow  from  the  blower  will 
start  only  when  the  temperature  of  the  battery  pack  rises  above  this 
set  value.  After  35  °C,  blower  will  come  in  to  deliver  5, 10  or  25  CFM 
of  cooling  air  per  module  (140,  280  or  700  CFM  for  a  battery  pack). 
The  battery  pack  delivers  power  for  176  min  (7.7  cycles),  69  min  (5.4 
cycles)  and  35  min  (3.5  cycles)  of  UDDS,  FIWFET  and  US06  driving 
cycles  respectively  before  reaching  the  cut  off  voltage  of  64.4  V.  The 
temperature  of  the  battery  pack  increased  constantly  during  the 
cycle  and  reaches  the  maximum  at  the  end  of  the  cycle.  The  heat 
generated  through  the  UDDS  and  FIWFET  cycle  is  lesser  and  hence 
can  be  effectively  dissipated  by  natural  convection.  Battery  module 
temperature  at  the  end  of  the  UDDS  and  HWFET  cycle  is  about  35  °C 
as  shown  in  Figs.  11  and  12.  As  shown  in  Fig.  13,  US06  cycle  has 
more  aggressive  driving  profile  as  compared  to  UDDS  and  HWFET 
cycles.  Hence,  more  energy  is  charged  and  discharged  to  and  from 
the  battery  pack  and  consequently  more  heat  is  generated.  The 
battery  pack  is  predicted  to  generate  4.93  MJ  of  energy  on  a  US06 
cycle  at  30  °C.  As  expected,  temperature  rise  of  the  battery  pack  is 
faster  in  US06  than  UDDS  and  HWFET  cycles.  The  blower  is  brought 
to  operation  at  744  s  when  the  average  temperature  of  the  battery 
pack  reaches  35  °C.  The  average  temperature  of  the  battery  module 
at  the  end  of  discharge  with  5  CFM  of  cooling  air  is  about  41  °C  and 
maximum  temperature  rise  at  the  end  of  the  cycle  is  about  11  °C.  In 
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Table  2 

Vehicle  and  cooling  system  specific  parameters. 


Parameter 

Value 

Vehicle  mass,  kg 

1828 

Frontal  area,  m2 

3.238 

Coefficient  of  drag,  Cd 

0.35 

Electric  motor 

75  1<W,  200  nm  max 

Battery  pack 

19.2  kWh 

Number  of  cell  per  module 

180 

Number  of  module 

28 

Amf  per  module,  m2 

0.00417 

Asm  per  module,  m2 

0.662 

Mass  flow  rate  of  cooling  air  for  battery  pack,  cfm 

140,  280,  700 

Surrounding  temperature,  °C 

30 

order  to  cool  the  battery  pack  further,  25  CFM  of  cooling  air  is 
blown  to  dissipate  the  heat  generated  and  average  temperature  of 
the  battery  module  is  brought  to  36.5  °C.  The  blower  is  operating 
for  about  22.65  min  of  the  total  duration  of  the  testing.  The  higher 
the  flow  rate  of  the  cooling  air,  the  slower  the  temperature  rise  of 
the  battery  module  but  parasitic  loss  and  power  consumption  will 
be  higher  for  higher  flow  rates.  Table  3  provides  the  total  energy 
used  for  propulsion  as  well  as  the  energy  recovered  through 
regenerative  braking  for  the  various  cycles.  UDDS  is  an  intensive 
start  and  stop  driving  cycle,  therefore  the  total  energy  used  and 
recovered  is  highest.  The  amount  of  energy  recovered  through  the 
HWFET  cycle  represents  the  highway  driving  condition  and  is  the 
lowest  among  the  three  driving  cycles.  The  selection  of  battery  cell 
to  build  a  battery  pack  is  an  essential  task  to  optimize  performance, 


Discharge  Capacity,  Ah 


Fig.  6.  Temperature  rise  of  the  battery  during  discharge  at  various  /t- rates. 


Fig.  8.  (a)  Current,  (b)  Voltage  during  the  first  6  min  cycle  of  the  SFUDS  profile. 


prolong  the  life  cycle  with  affordable  cost.  Smaller  cells  are  favor¬ 
able  for  thermal  management  and  cost  effective  but  high  number  of 
electrical  connections  may  result  in  higher  chances  of  failure  and 
energy  lost  especially  due  to  contact  resistance.  Although  larger 
cells  (prismatic  and  pouch  cell)  offer  lower  weight  and  volume 
efficiency,  cost  is  higher,  quality  of  the  cells  is  not  guaranteed  and 
difficult  in  terms  of  thermal  management.  Therefore,  a  careful 
consideration  is  required  to  select  the  type  of  cell  to  build  a  battery 
pack  for  a  certain  application. 

4.  Conclusions 


Empirical  equation  coupled  with  lumped  thermal  model  is  used 
to  predict  the  thermal  performance  of  the  LFP  cell  under  constant 
current  discharging  and  dynamic  charging  and  discharging  cycle. 


Time,  min 

Fig.  9.  Average  heat  generation  of  the  cell  during  600  min  of  the  SFUDS  cycle  predicted 
by  the  model. 
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Time,  min 


Time,  min 

Fig.  10.  Comparison  of  measured  and  simulated  SFUDS  profile  of  a  Li-ion  cell,  (a) 
Voltage,  (b)  Temperature. 


Fig.  11.  Thermal  response  of  battery  pack  to  UDDS  cycle  at  30  °C. 


Fig.  12.  Thermal  response  of  battery  pack  to  HWFET  cycle  at  30  °C. 


Time,  s 

Fig.  13.  Thermal  response  of  battery  pack  to  US06  cycle  at  30  °C. 


Table  3 

Energy  distribution. 


Driving  cycle 

Discharging,  MJ 

Charging,  MJ 

Heat  generation,  MJ 

UDDS 

81.91 

9.09 

0.108 

HWFET 

73.06 

2.17 

0.116 

US06 

58.99 

5.71 

0.176 

Good  agreement  between  the  simulated  results  and  experimental 
data  was  obtained.  The  heat  generated  from  the  cell  increases  with 
increasing  /rrates.  The  validated  battery  model  was  used  to 
investigate  the  evolution  of  the  battery  pack  temperature  of  a 
converted  EV  with  different  cooling  air  flow  rate  under  UDDS, 
HWFET  and  US06  driving  cycles.  The  heat  generated  from  the 
battery  module  is  highest  for  aggressive  US06  driving  cycle  and 
lowest  for  the  UDDS  driving  cycle.  In  a  less  aggressive  driving 
condition  for  UDDS  and  HWFET,  natural  convection  is  sufficient  to 
maintain  the  cell  temperature  at  35  °C.  On  the  other  hand,  US06 
cycle  requires  a  forced  convection  cooling  (19.0  W  m  2  K-1)  ach¬ 
ieved  by  blowing  700  CFM  of  cooling  air.  Hence,  a  well  designed 
thermal  management  system  is  needed  for  the  EV  battery  pack 
especially  under  aggressive  driving  conditions  to  ensure  safe  and 
reliable  operation  of  the  battery  pack. 
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